Revival of quantum correlations without system-environment back-action 
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Revivals of quantum correlations have often been explained in terms of back-action on quantum 
systems by their quantum environment(s). Here we consider a system of two independently evolving 
qubits, each locally interacting with a classical random external field. The environments of the 
qubits are also independent, and there is no back-action on the qubits. Nevertheless, entanglement, 
quantum discord and classical correlations between the two qubits may revive in this model. We 
explain the revivals in terms of correlations in a classical-quantum state of the environments and 
the qubits. Although classical states cannot store entanglement on their own, they can play a role 
in storing and reviving entanglement. It is important to know how the absence of back-action, 
or modelling an environment as classical, affects the kind of system time evolutions one is able to 
describe. We find a class of global time evolutions where back-action is absent and for which there is 
no loss of generality in modelling the environment as classical. Finally, we show that the revivals can 
be connected with the increase of a parameter used to quantify non-Markovianity of the singie-qubit 
dynamics. 

PACS numbers: 03.67.Mn, 03.65.Yz, 03.65.Ud 



I. INTRODUCTION 

Understanding the dynamics of correlations in quan- 
tum systems is essential for quantum information and 
quantum computation 0, and this has been thor- 
oughly investigated for quantum systems in quantum en- 
vironments. In this paper, we look at the behavior of 
correlations for the case when the environment is mod- 
elled as classical, and when the environment is unaffected 
by the system dynamics. This implies that there can be 
no back-action on the system via the environment. 

Sometimes only parts of the total correlations may be 
relevant. In particular, one can characterize, via appro- 
priate quantifiers, the quantum and classical parts of the 
correlations [H, Hj]. The quantum part of the correlations, 
quantified by quantum discord, includes not only entan- 
glement but also correlations that may occur in separa- 
ble states. For bipartite open quantum systems in quan- 
tum environments, entanglement may display phenom- 
ena such as sudden death, revivals and trapping 
Revivals of entanglement, after finite time periods when 
it completely disappears, can be expected when either 
direct two-qubit interactions @, Q or indirect effective 
interaction are present, as in the case of two qubits in a 
common quantum reservoir jToj . This is because interac- 
tions among quantum systems can create or destroy cor- 
relations, in particular, quantum correlations. Revivals 
may also occur for nonintcracting qubits in independent 
quantum environments and have then been related to 
the non-Markovian nature of the environments @, [Til ]. 
Recently, attention has been moved to the dynamics of 
correlations other than entanglement in the presence of 
either Markovian [I2J and non-Markovian [lj, [l4{ quan- 
tum environments. A peculiar aspect of the dynamics 



of these correlations is that even when total correlations 
behave in a regular way, both classical correlations and 
quantum discord may remain frozen for finite time inter- 
vals [H. 

In the case of initially entangled nonintcracting qubits 
in independent non-Markovian quantum reservoirs, en- 
tanglement revivals have been explained in terms of 
transfer of correlations back and forth from the two-qubit 
system to the various parts of the total system. This is 
due to the back-action via the environment on the sys- 
tem, which creates correlations between qubits and envi- 
ronments and between the environments themselves. In 
this paper we will refer to this kind of correlations as 
correlations induced by back-action. In particular, corre- 
lations may build u p b etween the two independent quan- 
tum reservoirs jlptilq and this phenomenon has been 
named sudden birth of entanglement in reservoirs. 

On the other hand, there are features that classical 
models of quantum systems fail to capture. Similarly, 
one would expect that if the quantum environment of 
a quantum system is modelled as classical, for example 
a classical light field coupled to a quantized atom, then 
there might be qualities of the dynamics that this model 
is unable to describe. Intuitively, one might expect that 
since a classical environment should not be able to store 
quantum correlations in the same way as a quantum en- 
vironment, the ability to describe revivals of quantum 
correlations might be affected. Similarly, if the environ- 
ment has no back-action on the system, then it would 
seem that this could also affect the correlation dynamics. 
Nevertheless, revivals of quantum correlations can occur 
for random classical telegraph noise [lj| . It is clearly im- 
portant to know exactly what kinds of behavior of corre- 
lation dynamics can be described when an environment is 
modelled as classical instead of quantum, or when back- 
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action on the system is not present. 

Existing work thus mostly concerns quantum systems 
in quantum environments, and has explained revivals in 
terms of back-action on the systems by their quantum 
environments, and in terms of quantum correlations in- 
volving the quantum environments. In this paper, we 
want to explain how and why revivals of correlations, in- 
cluding quantum correlations, can generally occur also if 
the environment is classical, and when no back-action is 
present. We illustrate our discussion with a model char- 
acterized by the absence of correlations induced by back- 
action, where the element of randomness is introduced 
in a very simple way. The example we will consider is a 
pair of independent qubits driven by single classical field 
modes with random phase, where revivals of correlations 
do occur both for entanglement, for the so-called quan- 
tum discord, and for classical correlations. 

We explain the revivals in terms of correlations in the 
classical-quantum state of the environments and qubits. 
We then discuss a class of global time evolutions for which 
back-action on the system is absent. This certainly oc- 
curs when the environment is unaffected by the system. 
We present a necessary condition for this to be the case. 
If the environment is initially in a classical state, this 
condition is also sufficient for it to be unaffected by the 
system. The state of a quantum environment may how- 
ever change. Interestingly, for this general class of time 
evolutions, modelling the environment as classical results 
in no loss of generality. In other words, any time evolu- 
tion of the system that could be obtained by coupling it 
to a quantum environment can also be obtained by mod- 
elling the environment as classical. Finally, we associate 
the revivals with the increase of a parameter that quan- 
tifies non-Markovianity of the reduced dynamics of the 
qubit. 

II. CORRELATION QUANTIFIERS 

Given a state p, we adopt the definitions in [2(| to 
quantify total correlations T(p), quantum discord D(p), 
entanglement E(p) and classical correlations C(p), 

T(p) = S{p\\ir p ) = S{tt p ) - S(p), 
D(p) ee S(p\\ Xp ) - S(x P ) - S(p), 
E(p) ee S(p\\a p ), 

C(p) ee S( Xp h Xp ) = S(ir Xp ) - S( Xp ), (1) 

where 5 , (p|jcr) ee — Tr(plogcr) — S(p) is the relative en- 
tropy, S(p) ee — Tr(plogp) is the von Neumann entropy; 
7Tp, X p an d o"p are respectively the product state, the 
classical state and the separable state closest to p, while 
tt Xp is the product state closest to xp- These states are 
such that they minimize the corresponding relative en- 
tropies. Note that entanglement E cannot be expressed, 
in general, as a difference of state entropies. The defini- 
tion of discord D(p) followed here and the original one 
Q in general do not coincide, but do so for Bell-diagonal 



states [2(| ■ We observe that the maximization procedure 
involved in the original definition of quantum discord has 
been analytically solved only for certain classes of quan- 
tum states, such as "X" states [2l[. Our analysis will be 
focused on the class of Bell-diagonal states, which is a 
subclass of "X" states, so that known analytical results 
[H H3 apply. 

For two independent quantum systems each interacting 
with its own local environment, so that neither direct nor 
mediated interactions between the two quantum systems 
exist, it is known that entanglement and total correla- 
tions cannot exceed their initial values [23|, [24| • Using the 
definition of T(p) in Eq. (p}, the definition of the prod- 
uct state closest to p and the "generalised H-theorem" 
for dynamical maps S(Ap\\Aa) < S l (p||cr) [25j], this prop- 
erty can be straightforwardly generalised to the case of a 
iV-partite system where each independent subsystem is 
locally affected by an arbitrary completely positive (CP) 
map. We note that in contrast to this, quantum discord 
can increase with respect to its initial value under local 
operations [2f|. 

III. MODEL 

We consider a pair of noninteracting qubits each locally 
coupled to a random external field, whose characteristics 
are unaffected by the qubit it is coupled to. This implies 
that back-action on the dynamics of the qubits is ab- 
sent. The aim is to give a clear example of how revivals 
of quantum correlations may occur even without back- 
action. Each environment is a classical field mode with 
amplitude fixed and equal for both qubits. The phase of 
each mode is not determined, and is equal either to zero 
or to 7r with probability p = 1/2 (the case when p\ ^ 1/2 
shall be treated in Sec. W\ . This model describes a spe- 
cial case of two qubits each subject to a phase noisy laser 
[27l l28j but where the phase can take only two values and 
with the diffusion coefficient in the master equation equal 
to zero. This model has in common with the one con- 
sidered in Ref. [l9| that the environments are unaffected 
by the system of interest. However, here the element of 
randomness is introduced in the phase which can assume 
only two possible values with assigned and fixed prob- 
abilities without switching between them. In Ref. [lj| 
qubits are subject to random telegraph noise, where in- 
stead it is the coupling that switches between two values 
during the evolution. 

For the case we are considering, the dynamical map 
for the single qubit S = A, B is of the random external 
fields type [8|, |25| and can be written as 

1 2 

As(i, 0)ps(0) = - ]T Uf (t) Ps (0)U- f (*). (2) 

i=l 

where Uf(t) = e~ lH i t / h is the time evolution operator 
with Hi = \hg(a + e^ 1< ^ i — er-e 1 ^) and the factor 1/2 arises 
from the equal field phase probabilities of the model 
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FIG. 1: (Color online) Sketch of the physical system. The two 
qubits are initially entangled. Each random dephaser (RD) is such 
that the field phase at the location of each qubit is either zero or ix 
with equal probability, p = 1/2, corresponding to either Eg or EjJ 
(S = A, B). 



(more in general, there is a probability pf associated to 
Uf ). Each Hamiltonian Hi is expressed in the rotating 
frame at the qubit-field resonant frequency uj. In the ba- 
sis {|1), |0)}, the time evolution operators Uf(t) have the 
matrix form 



uf(t) 



cos(gt) e sin(gt) 
-e 1 ^ sin(gt) cos(gt) 



(3) 



where i = 1, 2 with <f>\ = and 4>2 = tt- 

This simple model is depicted in in Fig. [1] where the 
two local environments are realized by two single clas- 
sical field modes of the same frequency and amplitude, 
but each passing across a random dephaser such that the 
phase of the field at the location of each qubit is either 
zero or 7T, with probability p = 1/2. The interaction 
between each qubit and its local field mode is assumed 
to be strong enough so that, for sufficiently long times, 
the dissipation effects of the vacuum radiation modes on 
the qubit dynamics can be neglected. The dynamics is 
cyclic, as evinced from the fact that the unitary matrix 
of Eq. ((3|) (for each i = 1,2) becomes the identity for 
times t = kw/g (k = 0, 1, 2, . . .). The same initial state is 
thus retrieved at these times. If one takes into account 
dissipation effects on the system, the dynamics would 
not be cyclic anymore, but it is expected that for dissi- 
pation weak enough the qualitative behavior of the dy- 
namics would remain the same until a characteristic time. 
This could be feasible by considering as qubits two atoms 
placed in separated zero temperature cavities. Each atom 
is subject to a resonant interaction with a laser with a 
random phase but out of resonance with respect to the 
characteristic cavity frequencies in order to slow down 
spontaneous emission. 

The global dynamical map A applied to an initial state 
p(0) of the bipartite system, p(t) = A(t, 0)p(0), is again 
of the random external fields type, that is, 



= I E U?{t)U?{t)p{Q)U?\W!\t)- (4) 



We now show that the map A moves inside the class 
of Bell-dia gon al states (or states with maximally mixed 
marginals j22[). We first introduce the notation \1±) = 



(|01) ± |10))/V2 for the one-excitation Bell states and 
|2±) = (|00)±|ll))/ v / 2for the two-excitation Bell states. 
By using Eqns. ([3]) and ((4]), we find that the map acts on 
single Bell states as 

A(t,0)\i±){i±\ = [1 - /(*)]l*±><i±l + f(t)\i' T ){i' T l (5) 

where = 1,2 with i ^ i' and f(t) = sin 2 (2gi)/2, so 
that mixtures of the states |1 + ) (|1-)) and |2_) (|2 + )) 
are mapped onto different mixtures of the same states. 
It immediately follows that the map A moves inside the 
class of Bell-diagonal states, with a generic initial Bcll- 
diagonal state written as 



P (o) = T / xi(o)\i s )(i l 



i,s 



|, (i=l,2;« = ±). (6) 



It is worth noting that the diagonal Bell states in 
Eq. © include the well-known Werner states [ll[ and 
that p(0) is entangled if the largest Xf (0) is greater than 
1/2 [20]. The action of the map (UJ) on p(0) thus gives 
another Bell-diagonal state p[t) with time-dependent co- 
efficients 

Xt(t) = Xt(0)[l-f(t)]+X^(0)f(t) (i^i')- (7) 

It can be shown that, for all t > 0, the largest Af (t) is 
always the same, while, depending on the initial values 
of the coefficients, there are switching times after which 
the second largest Xf(t) changes. Moreover, since the 
function f(t) is periodic, there will be repeated switch- 
ing times during the dynamics. The restriction of the 
dynamical map of Eq. Q to the class of Bell-diagonal 
states and the order between the Af (t) is crucial to con- 
struct the classical, separable and product states closest 
to p(t) and thus to analytically calculate the correlation 
quantifiers in Eq. ([T]) [2fJ. In particular, the closest prod- 
uct states do not evolve (7r t = A7Tn = ttq) and arc all equal 
to the normalized 4x4 identity 1/4, while the classical 
state closest to p(t) is equal to the evolved classical state 
closest to p(0), that is x P (t) = Xa p (o) = A X P (o)- These 
properties permit us, by using the generalised H-theorem, 
to show that 



D[p(t)] = S(p(t)\\ Xm 
C[p(t)}=S( Xpi t)\h 



< S(p(0)\\ Xm ) = D[p(0)}, (8) 
, pW )<S(p(0)h P (.o)) = C[p(0)}. 



Therefore, under the map of Eq. ((4} not only entangle- 
ment but also discord and classical correlations cannot 
exceed their initial values during the dynamics. 



IV. DYNAMICS OF CORRELATIONS FOR 
EQUAL FIELD PHASE PROBABILITIES 

We are now able to analyze the evolution of the rele- 
vant quantifiers of correlations under the map of Eq. (U) 
for different choices of p(0). One finds a rich variety of 
behaviors of the quantifiers. In particular, Fig. [2] shows 
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FIG. 2: (Color online) Dynamics of discord D (blue solid line), 
classical correlations C (red dashed line) and total correlations T 
(green dotted line) for an initial Bell-diagonal state with A^(0) = 
0.9, A~(0) = 0.1 and \f(0) = Aj (0) = 0. Total correlations oscil- 
late regularly, while transitions between decay and rise of quantum 
and classical correlations occur periodically. 





\ 

1 


TV 




if 




\ 

\ 

\ 

\ 

\ 

\\ 




\ 

\ 

\ 

\ 

\ 


1 1 
1 1 



)i v . .«y ^ >y i 

f, =7r/4 nil t 2 =3n/4 n 



gt 

FIG. 3: (Color online) Dynamics of entanglement E (orange solid 
line) and discord D (blue dashed line) for an initial Bell-diagonal 
state with the same coefficients as in Fig. [2] D vanishes only at 
given times t„ = (2n — l)n/(4g) (n = 1,2,...), while E disap- 
pears for finite time intervals. The quantifier of non-Markovianity 
X E {t)/ 10 (green dotted line) is also plotted, showing its increase in 
correspondence with the revival of correlations. 



the time evolutions of T, D and C starting from a Bell- 
diagonal state consisting of only the two Bell states |1±), 
with \f(0) = 0.9, Ar(0) = 0.1. This is a case for which 
switching times for the time-dependent coefficients (t) 
and Xi (t) exist. In this case the dynamics of correlations 
is such that, while total correlations oscillate, either clas- 
sical correlations or discord are frozen during alternat- 
ing finite time periods. When discord is frozen, classical 
correlations are oscillating, and vice versa. This kind of 
behavior has already been found in the case of two qubits 
locally subject to non-dissipative channels [l5| • This may 
thus be taken to signify that such peculiar behavior of 
both quantum and classical correlations is a general re- 
sult not connected to specific systems and interactions. 

Fig. |3] displays the dynamics of entanglement E start- 
ing from the same initial Bell-diagonal state as before. 
Again, collapses and revivals of entanglement occur, as 
is known to happen for quantum non-Markovian envi- 
ronments such as amplitude damping channels @. In 
our case, the occurrence of revivals is clearly connected 
to the cyclical nature of the total dynamics. We note that 
in presence of small dissipation, the same qualitative be- 
havior of correlations is expected, that is, partial revivals 
should still occur when the dissipation effects are weak 
enough. We also note that the occurrence of entangle- 
ment revivals is not in conflict with the fact that the two 
qubits arc subject to local independent environments. In- 
deed, while entanglement must be constant if the local 
actions are unitary, in the general case for non-unitary 
dynamics (as in our model) the only constraint (already 
mentioned above) is that the entanglement cannot ex- 
ceed its initial value during the dynamics. This is not 
in conflict with the fact that, for example, entanglement 
goes to zero and then revives up to the initial value. 

It should be pointed out that the above revivals of 
quantum correlations, as characterized by entanglement 



and discord, occur without transfer of correlations in- 
duced by back-action from system to environment. The 
absence of back-action implies that the initial two-qubit 
correlations cannot be "temporarily stored" as correla- 
tions entirely within the environments. This differs from 
a common point of view where revivals are interpreted in 
terms of transfer of correlations induced by back-action 
which correlate the independent non-Markovian environ- 
ments, which in turn again re-correlate the qubits fltl - 
[l8j . In our example this mechanism is absent. 



V. DYNAMICS OF CORRELATIONS FOR 
DIFFERENT FIELD PHASE PROBABILITIES 



The above results are obtained in the case when the 
phase of each field is zero or n with equal probability 
distributions for both the subsystems (pf = p§ = 1/2, 
S = A,B). Extending the map of Eq. ((4J to the case 
when pf = 1 — p§ 7^ 1/2 it becomes 



AMMO) = £ ptpfU?{t)U?{t)p{Q)U?\t)Uf\t). 

(9) 

This map, applied to an arbitrary Bell-diagonal state of 
Eq. ([5]), in the standard computational basis B = {|1) = 
|11),|2) = |10),|3) = |01),|4) = 1 00)} gives the state at 
time t 



m 



fa(t) 

b(t) 



b(t) 

I - a(t) 
c{t) e(t) 
\d(t) -c(t) 



c(t) d(t) \ 

e(t) -c(t) 

-a(t) -b(t) 

-b(t) a(t) J 



(10) 
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where 

a(t) = 

b(t) = 

c(t) = 

d(t) = 

e(t) = 



[1 - A+(0) - Ar(0)]/2 - [-F(A) - G(X)I(p)]f(t), 

[F(X) - pfG(X) - pfL(X)] sin(4^)/4, 

[F(A) - pfG(X) - pfL(X)} sm(4fft)/4, 

[F(X)-I(p)L(X)]f(t)-Xz/2, 

d(t) + [l-2Ar(0)-A+(0)]/2, (11) 



with F(X) = 1 - 2A+(0) - Ai (0) - A+(0)/2, G(A) = 
1 - 2A+(0) - A+(0), L(A) = 1 - 2[A+(0) + A^O)], 
lip) = P?+P?-2PiP? and /(*) = sin 2 (2.gi)/2, as before. 
The density matrix p(t) above is not Bell-diagonal in the 
basis B [22j , but it is straightforward to show that it has 
maximally mixed marginals, Tr^/5(t) = Tr^p(t) = 1/2, 
just like Bell-diagonal states. In fact, it is known that 
any two-qubit state with maximally mixed marginals al- 
ways can be reduced, up to local unitary equivalence, to 
a state which is Bell-diagonal in some other basis [22l. r29l] . 
Since local unitary operations do not affect the correla- 
tion quantifiers Q , we can use the expressions of Eq. ([T]) 
to study the dynamics in terms of the eigenvalues of 
p(t), analogously to what was done above for the case 
pf = 1/2. Notice that, as a consequence, classical corre- 
lations can be obtained by C = T — D, a relation valid 
only for Bell-diagonal states (20j . 

For the sake of simplicity we limit the present analy- 
sis to the case when pf = pf = pi, because the other 
cases Pi ^ pf only result in quantitative but not qual- 
itative differences in the time behavior of the correla- 
tion quantifiers. We are able to plot total correlations 
T(p(t)), discord D(p(t)) and entanglement E(p(t)) as 
functions of gt by choosing suitable values of the ini- 
tial coefficients A| (0) and diagonalizing the state p(t) for 
some values of p\ at different values of t. In Fig. @] we 
in particular display the plots for the same initial Bell- 
diagonal state as was chosen for the case pf = 1/2 (that 
is, Xf (0) = 1 — Aj7(0) = 0.9 and A^(0) = 0). It is 
seen that the amplitudes of oscillation for all the correla- 
tion quantifiers tend to decrease when p\ moves from 1/2 
(the behavior is symmetric with respect to 1/2), with the 
quantifiers remaining closer to their initial values. In par- 
ticular, quantum discord and entanglement do not vanish 
anymore. This behavior is expected, because changing p\ 
cannot affect the presence or the frequency of oscillations, 
as is clear from the cyclic behavior of the dynamics due 
to Eqs. §5$ and ©■ Moreover, when the limiting case 
is considered where each local field becomes fully deter- 
ministic (pf = 0,1), the two-qubit evolution is a tensor 
product of local unitary evolutions on each qubit (see 
Eq. ©), both when pf = pf (fields in phase) and when 
Pi = 1 — pf (fields with opposite phase). In this case all 
correlation quantifiers are, as known Q , constant in time 
and equal to their initial value: the only possible behav- 
ior for correlation quantifiers when p\ goes to (or 1) 
is thus a reduction of the oscillation amplitude towards 
their initial value. 

These findings show how the most relevant case for 
our purposes is the case when pf = = 1/2, when both 
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FIG. 4: (Color online) Dynamics of correlation quantifiers T 
(dashed lines, panel a), D (solid lines, panel a) and E (panel b) 
for an initial Bell-diagonal state with the same coefficients as in 
Fig. [2] and for values of pi = 1 — p2 from bottom to top equal to 
0.2 (orange line), 0.08 (red line), 0.025 (blue line). Discord D is 
almost constant for p± = 0.025. 



entanglement and discord vanish and revivals occur. 



VI. CORRELATIONS IN A 
CLASSICAL-QUANTUM STATE 

We have shown that quantum correlations can revive 
also for a classical environment, and without back-action. 
To look more closely at this, let us consider a so-called 
classical-quantum state of the two qubits labeled A and 
B, and the classical environments labeled a and b. This 
state is initially given by 



2 2 

PQE (0) = p% B (0) ®$>H)H<i|®£ Pi \j) b E b E (J 1,(12) 

i=l 7=1 



where Pq B is the quantum state of the two qubits, and 
|1)^ and |2)^ are two "classical" basis states that envi- 
ronment a can take, and |1)^ and \2) E are classical basis 
states of environment b. The classical environments can 
only be in separable mixtures of their basis states; in our 
case, pi — P2 = 1/2. The single-qubit map on qubit A is 
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now equivalent to applying a unitary transform 
2 

u Aa {t) = Y,uf{t)®\i) a M (is) 

i=l 

to qubit A and its environment a, where U A is defined 
in Eq. ([2J for S = A, and tracing over the environment. 
The environment acts as "control system" for what uni- 
tary operation is applied to the qubit. The map U Bh (t) 
for qubit B is constructed analogously. During time evo- 
lution according to U Aa and U Bb , the states of the clas- 
sical environments do not change. The qubits obviously 
do not affect their environments, and back-action by the 
environments on the qubits is clearly not present. The 
two-qubit map in Eq. Q can be obtained by applying 
U Aa (t) ® U Bb (t) to both qubits and environments, and 
tracing over the latter. 

We can now quantify the correlations in the classical- 
quantum system exactly as for a quantum-only system, 
using the definitions in Eq. ([l}. Since the subsystems 
Aa and Bb both evolve under local unitary operations, 
the quantum and classical correlations between Aa and 
Bb are obviously constant. If we trace over the environ- 
ments, then the correlations that remain between A and 
B were found to oscillate. If we trace over the qubits A 
and B, then the environments a and b remain completely 
uncorrelated throughout the evolution. 

If we look at the qubit-environment state Aa, "entan- 
glement" arises neither within this state, nor in Bb. In- 
terestingly, "quantum discord" may arise in the classical- 
quantum state Aa, and in Bb, however, one easily realizes 
that the presence of "quantum discord" is not necessary 
for the revival of quantum correlations between the two 
qubits (32l | . The classical environments cannot store any 
quantum correlations on their own, and they do not be- 
come entangled with their respective quantum systems. 
But through the classical record they keep memory of 
what unitary operation has been applied to the qubits, 
they can obviously nevertheless play a vital part in re- 
viving the quantum correlations, in particular the en- 
tanglement, that was initially present in the quantum 
systems A and B. This remains the case even without 
back-action. Although the environments can hold infor- 
mation about the quantum systems, it may be somewhat 
misleading to talk about "information flow from the sys- 
tcm(s) to the environment (s) " , or "information backflow 
from the environment (s) to the system(s)". This is be- 
cause the information an environment holds about a sys- 
tem is due to what action the environment has on the 
system, and not the other way around. 

After this explanation, it may seem trivial that classi- 
cal environments can store quantum correlations in the 
sense that the state Aa — Bb can be entangled, but the 
state A — B when the classical environments are traced 
out is less entangled or even separable. This fact has 
nevertheless been overlooked when explaining the mech- 
anisms for revivals of quantum correlations. Although 
in our case the environments are classical, quantum cor- 



relations in the system can obviously revive when there 
is no back-action also for quantum environments. The 
considerations above can also easily be generalized to ar- 
bitrary number of dimensions for both systems and en- 
vironments. 



VII. A CLASS OF GLOBAL TIME 
EVOLUTIONS GIVING NO BACK-ACTION 

It is now clear that classical and quantum correla- 
tions may revive even when environments are modelled 
as classical, and for no back-action, and why this hap- 
pens. There is no back-action when the system does 
not affect itself via the environment. This must be the 
case when the environment dynamics is unaffected by the 
system. We will now discuss a class of global (system- 
environment) evolutions where back-action is absent. We 
show that for this class of time evolutions, any reduced 
dynamics of the system that can be obtained with a quan- 
tum environment, can also be obtained with a classical 
environment. 

We consider time evolutions (in the finite-dimensional 
case) described by completely positive (CP) maps pj. 
Such time evolution of a quantum system S can be de- 
scribed as a unitary transform acting on S and some en- 
vironment E, followed by tracing over E. This will yield 
a CP map for the time evolution of S, 

Ps{t) = Tr E [U(t)p s (0)®p E (0)UHt)} 

= ^A k (t)p Q (0)Al(t). (14) 

k 

Here Afc(t) are the Kraus operators for the CP map, sat- 
isfying ^2 k A' k (t)Ak{t) = Is, and pe{0) is some initial 
state of the environment. (Given a CP map, we can con- 
struct U(t) so that p E (0) = |0)bb(0|.) 

The unitary transformation U(t) will in general also 
affect the state of the environment so that not all com- 
pletely positive maps on a quantum system can be de- 
scribed trough unitary evolution of the system plus an 
unchanged environment. Denote the basis states for S 
and E by \i)s and \j)e, respectively. We also assume 
that a classical environment could be prepared in any 
mixture of its basis states, and a quantum environment 
in any linear superposition of these states. We moreover 
let the environment undergo an assigned evolution UE{t) 
which is unaffected by the system, meaning that its evo- 
lution is just pE(t) = UE{t)pE{0)U E (t). One realizes 
that the condition 

s(i\E(j\U E (t)U(t)\l) E \k) s = unless j = I (15) 

is necessary for the state of the environment to be unaf- 
fected by the system. It is also possible to show that this 
condition is sufficient for the environment not to be af- 
fected by the system provided that it is initially prepared 
in a classical mixture of the chosen basis states. This 
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can be shown by looking at pE(t) and, using the condi- 
tion in Eq. (|15p . requiring that it satisfies the equality 
p E (t) = U E {t)p E (0)U E (t). If Eq. flTSJ) is satisfied, then 
U(t) must be of the form 



u(t) = Y, u ^ t )® u E(m)E E {]i 



(16) 



where U?(t) = E(j\U E (t)U (t)\j) e are arbitrary unitary 
transforms on the system. The time evolution of the 
quantum system can therefore be described by the ran- 
dom application of unitary transformations, 



Ps(t) 



Y,pjuS(t)p S {w 3 



(17) 



where the time-independent probabilities pj are deter- 
mined by the initial state of the environment. This state 
can be chosen as a quantum state \iP)e = J2j y/Pj\j)E, 
but the "classical" state pe = J2j Pj\j)EE(j\ will yield 
exactly the same CP map for the system. Importantly, 
therefore, if the global time evolution has the form of 
Eq. (p~6|) . there is no loss of generality in modelling the 
environment as classical. It is worth noticing that the 
map acting on the system given in Eq. (|17[) is exactly of 
the "random external fields" type. The map considered 
in our physical model in Sec. IIIII is of this form. 

We remark that if the environment is in a quantum 
superposition state of the relevant basis states, then the 
reduced state of the environment can be affected by the 
system, since in this case the system and environment 
may become entangled. More specifically, diagonal ele- 
ments of the environment reduced density matrix evolve 
unaffected by the system, according to L/s(t) only. Non- 
diagonal elements, however, including their absolute val- 
ues, can be affected by the system dynamics. Neverthe- 
less, this behavior in the environment reduced density 
matrix has no effect on the reduced dynamics of the sys- 
tem. 



VIII. REVIVALS AND NON-MARKOVIANITY 

It also looks appealing to interpret the revival of cor- 
relations by referring to intrinsic characteristics of two- 
qubit system evolution itself. Revival of correlations can 
be connected to the degree of non-Markovianity of the 
map, which effectively is the memory that the state of 
the system has of the correlations it owned in the past. 
In our model of the two-qubit dynamics, negative rates in 
the corresponding single-qubit master equation appear, 
as occurs e.g. for a phase noisy laser [28| . of which our 
model is a special case. Negative rates typically indicate 
non-Markovianity of the map [3(| ■ Both the single-qubit 
and the two-qubit time evolution in our model is thus 
non-Markovian. In particular, in Fig. [3] we see that the 
state p(ti) at ti = n /4g has both zero entanglement and 
zero discord. Nevertheless, both entanglement and dis- 
cord subsequently revive. If we take the intermediate 



state p(t\) as initial state, then since correlations cannot 
exceed their initial values as shown by Eq. ©, this state 
remains non-entangled and has zero discord. Therefore, 
A(Mi)p(ii) + A(t,t Q )p(t ) where p{t x ) = A(t 1 ,t Q )p(t ). 
Our map does not satisfy the Markovian composition 
law A(tAo) = A(t, ii)A(ti, to) for completely positive 



maps 25[, where A(t2,ii) has been defined in Eq. (QJ. 



This means that the time evolution from the intermedi- 
ate state p(ti) onwards is not a CP evolution. 

One can render the association between revival of cor- 
relations and the non-Markovianity of the time evolution 
quantitative as follows. We generalize the quantifier of 
non-Markovianity introduced in Ref. [3l[ , which is based 
on the increase of entanglement between the system and 
an isolated ancilla, to the time dependent form 



T E {t) 



to 



dE[p AB {t')} 



df 



dt'-AE(t), (18) 



where E is defined in Eq. (P) and AE(t) = E(t) - 
E(to), with t < t. Eq. (fTHj) can be applied to the 
single-qubit map ([2]) in our model by initially prepar- 
ing qubit B (system) and another qubit A (ancilla) in 
a maximally entangled state, for example pab(0) = 
|2 + )(2 + |. The evolved state is then given by pab {t) = 
cos 2 (^)|2+)(2 + | +sin 2 (gt)|l_)(l_| while the closest sep- 
arable state (T pAB (t) = ^[|2 + )(2 + | + |1_)(1_|] is time in- 
dependent. We find that Z E (t) monotonically increases 
in the same time regions when revivals of entanglement 
and quantum discord occur, as shown in Fig. |3l As a 
consequence of this behavior, X E (t) would increase to in- 
finity in the ideal case of zero dissipation. In a more 
realistic scenario, the periodic behavior of various quan- 
tifiers will be limited to a small number of periods and 
I E (t) will increase until dissipation destroys the revival 
of correlations. 

The entanglement in Eq. (TT5)) is that between a single 
qubit subject to some single-qubit time evolution, and 
an ancilla qubit which does not evolve. One can also 
define a corresponding quantifier for which the ancilla 
qubit is also subject the same single-qubit time evolution 
as the first qubit. The entanglement in this quantifier 
would then correspond exactly to that between the two 
qubits in our model, as plotted in Fig. [3] This makes 
the connection between the revival of two-qubit system 
correlations and the non-Markovianity. 



IX. CONCLUSIONS 

We have considered a simple model where random 
external classical fields act locally on two independent 
qubits. We found that collapses and revivals of entan- 
glement between the qubits, as well as decreases and 
increases of quantum discord and classical correlations, 
occur even when the environment is unaffected by the 
system qubits. This is in contrast to the common inter- 
pretation that revival of correlations is caused by trans- 
fer of back-action-induced correlations between system 
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and environment [161 1 1 81 ] . In terms of correlations in a 
classical-quantum state, we have explained why quantum 
correlations can revive also for classical independent en- 
vironments on which back-action is not present. We have 
shown that the increase of correlations in our system can 
be connected with the increase of a parameter used to 
quantify non-Markovianity of the dynamics of the indi- 
vidual qubits. If randomness would be introduced in the 
amplitude instead of in the phase, as considered in this 
paper, one may expect that revivals of correlations should 
also occur. 

Furthermore, we have discussed a class of global time 
evolutions for which back-action is absent, even if the 
reduced dynamics of the environment can be affected by 
the system. We have shown that for this class of time evo- 
lutions, any system dynamics obtained with a quantum 
environment can also be obtained by modelling the envi- 
ronment as classical. This global time evolution results in 
a reduced dynamics for the system precisely of the form 
that arises from random external fields. Suppose that 



a quantum system plus its environment undergo unitary 
evolution, so that there is no additional environment. It 
is natural to think of a classical environment as a "macro- 
scopic" system that should not be significantly affected 
by the "microscopic" quantum system it interacts with. 
Interestingly, what we found means that conversely, if a 
quantum system either does not affect its environment 
or influences the latter in a way that does not result in 
back-action, then there is no loss of generality in mod- 
elling the environment as classical, as far as the possible 
time evolutions of the system are concerned. 
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